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Abstract. Random optical-pulse polarization switching along an active optical medium 
in the A-configuration with spatially disordered occupation numbers of its lower energy sub- 
level pair is described using the idealized integrable Maxwell-Bloch model. Analytical results 
describing the light polarization-switching statistics for the single self-induced transparency 
pulse are compared with statistics obtained from direct Monte-Carlo numerical simulations. 



The model of light interacting with a material sample composed of three-level active atoms 
has made possible the descriptions of several nontrivial optical phenomena, including lasing 
without inversion [TT], slow light [12], and electric-field polarization of solitons in self-induced 
transparency [13j. Its simplest version including a non-degenerate upper and two degenerate 
lower working atomic levels — the A configuration — is completely integrable when the 
pulse width is much shorter than the medium relaxation times [H]. It describes a new type 
of a self-induced transparency pulse, which may be a solitary wave only asymptotically, but 
in general switches into one of the two purely two-level transitions between one of the lower 
levels and the upper level. These transitions correspond to circularly polarized light, and 
the direction of the switching is determined by the population sizes of the degenerate lower 
levels [ini [16]. Thus, for spatially disordered populations, random polarization switching 
takes place as a light soliton travels along the material sample. The integrability of the A 
configuration furnishes a unique opportunity to study the mechanism responsible for this 
random switching and its statistical properties exactly in the framework of a sufficiently 
idealized model, which otherwise would be impossible because of strong nonlinearity. In this 
letter, we both discuss the analytical results [T7] on this switching and compare them with 
the results of numerical simulations. 

Resonant propagation of ultra-short, monochromatic, elliptically polarized light pulses 
through a two-level, active medium with a doubly degenerate ground level (A-configuration) 
is described by the quasi-classical Maxwell-Bloch system [181 El 1131 IIS] 

/oo 
-oo 

(lb) dtp+ - 2i\p+ = [E+ {U-n+)- E^p*] /2, 

(Ic) dtp- - 2t\p_ = \E_{U - n_) - E+p] /2, 
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(Id) 9,/i=[E+V+i5;_p+*]/2, 

(le) dtU = - + E+*p+ + E^pJ + EJpJ\ /2, 

(If) dtn^ = [E±p±* + /2. 



Here, E±{x,t) are the envelopes of the electric field and p±{x,t,\) and p{x,t,X) of the 
medium-polarization, n±{x,t,X) and Af{x,t,X) the population densities of the ground and 
excited levels, respectively, A the frequency detuning, and g{X) > 0, with g{X) dX = 1, 
the spectral-line shape. The and "— " transitions interact with the left- and right- 
circularly polarized pulse components, while p is due to the two-photon transition between 
the ground levels. The purely two-level "+" and "— " transitions are invariant and involve 
only circularly-polarized light. A time-conserved quantity of Eqs. ([T]) is N' + + n_ = 1, 
where unit normalization is chosen. 

The approximations made in Eqs. ([T]) are (i) the pulse-width is much longer than the light 
oscillation period (slowly- varying envelope approximation) and much shorter than the relax- 
ation time-scales in the medium, and (ii) unidirectional propagation. The latter holds pro- 
vided the interaction time of counter-propagating pulses is much shorter than the nonlinear- 
response time of the medium. Equations ([1]) are dimensionless, e.g., the speed of light is 
c = 1. 

If the spectral width of the pump pulse priming the ground states of the medium is much 
broader than the width of g{X), i.e., the initial populations can be considered homogeneous 
within the width of g{X), we find the two components of the soliton solution [151 [ISl [I7| 

2(3{t -x) + TX 

(2) + I In + i In cosh (2rA{x) + In ^) " , 
where 

G±{x) = v/[l±tanh(2rA(x) +ln|rf+|/|d_|)] /2 
are their amplitudes and 

6±(x, t) = 2'y(t — x) + (7[x ± A{x)] — a.Tgd± 
their phases, which in turn depend on the soliton parameters 7 and /3. Here, 

(3) A{x) = j; aiOd^, 

is the cumulative initial population difference a{x) along the medium sample up to any given 
position X, which satisfies the asymptotic condition 

(4) lim n±{x, t, A) = [1 ± a{x)] /2 > 0. 

t— >— 00 

The rest of the material variables are known to vanish as t — t- —00 for this solution [16], so 
that only the two degenerate lower levels are populated initially. Putting the initial time at 
—00 is justified because, in gases, the lifetime of the system ranges from 10~^ to 10~^ seconds, 
while the typical pulse- width is 10~^ seconds or shorter [19]. The real- valued coefficients a 
and T are given by 



E±{x,t) = 42/3G±(x)e*®±(^'*) sech 
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with /3 > 0. Equation ([2]) shows that the maximal amphtude of each sohton component 
equals 4/3 and its temporal width equals 1/(2/3). The constants d± give the soliton phase 
and position. Note that, since r < 0, the amplitude decreases and G-{x) increases 

with increasing A{x), and vice versa with decreasing A{x), which is the polarization-switching 
effect of [ISlIinj. 

The light-pulse polarization can be described in terms of the polarization ellipse, which 
is characterized by the orientation and ellipticity angles, ip ^ind t], with — 7r/4 < rj < 7r/4. 
These can be found from the formulas 

tan2tp = i{E+E*_ - E_El)/{E+E*_ + E^E^) 

and 

sin2r7 = (|E+|2-|E_|2)/(|E^|2 + |i5;_|2), 
which, for the sohton (j2]) give [16] 

^ = -(yA{x) + arg /2, 
sin2r/ = tanh [2tA{x) + In / . 

Note that these two angles are time-independent. 

If the initial population difference a{x) in the medium is random and spatially statistically 
homogeneous, we can approximate it as white noise 

(6) (a(a;)) = 6, i\a{x) - h]{a{x') - h]) = a^5{x - x'), 

where (■) denotes ensemble averaging over all possible realizations of and is the 

Dirac Delta function. This approximation is consistent provided the pulse-carrier frequency 
Aq, the correlation length of the soliton width 1//3, and the observation location x 
along the sample satisfy the inequalities \q ^ ^ Ij P <^ x. The first is related to the 
slowly- varying envelope approximation (mentioned above), the second to the unidirection- 
ality assumption, and the last to the white-noise assumption. In this case, the cumulative 
integral A{x) in Eq. ([3]) can be approximated as A{x) ~ aW{x) + hx, where W{x) is the 
standard Wiener process. Note that the parameter a in Eq. ([6]), the correlation length Lc, 
and the variance cr^ of a true initial population difference a{x) are related by a = \/2Lc a a- 

In an experiment, Lc would be approximately the same as the coherence length ic ~ 
A^/AAp of the pump light used to prepare the optical medium, where Ap is the average 
wavelength of the pump light and AAp the characteristic width of the light-source spectral 
line. If the pump was a Ti-sapphire laser, Ap ~ 800 nm and AAp ~ 5 nm [20], so lc ~ 
0.1 mm ^ Aq (~ 600 nm for sodium vapor) and a several-centimeters long experimental 
device would be sufficiently long to capture the desired statistical effects. 

The orientation angle il){x) behaves like a Brownian motion with drift —ah and diffusion 
coefficient ^cr^a^and its probability density function (PDF), p.^{x; s), at any x is Gaussian in s 
with mean {ip{x)) = —(rbx+^ axg{d*id+) and variance cr|,(x) = a'^a^x. Note that the value of 
ip{x) is fixed at ip = arg((ilci+)/2 when cr = 0. The PDF p^{x; s) shows excellent agreement 
with numerical simulations in Fig. [H The Lorentzian spectral-line shape g{X) = e/7r(A^ + £^) 
was used in determining p^{x; s) and all subsequent PDFs. Comparisons in all figures are 
made with e = corresponding to the Dirac Delta function, i.e., the limit of an infinitely 
sharp spectral line. 
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Figure 1. PDF p^{x;s), with (3 = 1/3, 7 = 1/3, e = 0, d+ = d_ = i, 
theoretical (black lines) and results from 1600 simulations (gray lines; green 
online). Left: b = 0,a = 0.75. Right: b = -0.75, a = 0.5. 
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Figure 2. PDF Pr^ix;s), with (3 = 1/3, 7 = 1/3, e = 0, d+ = d_ = i, 
theoretical (black lines) and results from 1600 simulations (gray lines; green 
online). Left: b = 0,a = 0.75. Right: b = -0.75, a = 0.5. 



The PDF for the ellipticity angle rj at any x equals 

Pri{x; s) = (1 /V271X a|r| cos 2s) exp{ — [tanh^"^(sin2s) — 2r6x — In |(i+|/|c?_|]^/8a^r^x} 

for —tt/A < s < 7r/4. For large x, Pyix; s) concentrates at one or both circular polarizations 
at s = ±7r/4 clS p^yj i^X] s ) ~ [6{s + 7r/4) + 6{s - 7r/4)]/2 if 6 = 0, ~ 6{s + if 6 > and 

~ 6{s — 7r/4) if 6 < 0. The PDF Pri{x; s) shown in Fig. [2] accurately describes the numerical 
simulations. 

To quantify the polarization switching statistics, we consider the periods during which 
the ellipticity angle 77 enters, exits, and stays in the vicinity of either circular polarization, 
7r/4 — K = < |?7| < vr/4, for some appropriately chosen small angular distance k. More 
precisely, we study the statistics of the following random distances: Xtra, over which the 
pulse polarization evolves from linear [r] = 0) to nearly-circular of either orientation {\r]\ = 
7]^); Xint, over which the pulse polarization evolves from linear to nearly-circular of either 
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Figure 3. Comparison of exit time statistics between analytical formulas and 
results from numerical simulations when (left) no bias, b = 0, a = 0.75, and 
(right) b = -0.75, a = 0.5. For both, (3 = 1/3, 7 = 1/3, e = 0, d+ = d. = i. 



orientation and back to linear; and Xdep, beyond which the pulse polarization remains forever 
circularly-polarized, \r]\ > rj^^, for all greater distances. 
When b = (a(x)) = we find [17j 

00 

(7a) px,Jx) = ^/2L/^ Yl + l)e-(^"+l)'^/2^ 

(7b) px,Jx) = ^/2L/^ 

where L = [tanh~^(cos2K)/2a|r|]^. These two distributions are depicted in Fig. El^left) and 
agree with the numerical simulations. Note that (Xtra) = L, cr\^^_^ = 2L^/3, but (Xint) = 00. 

Another statistic, when 6 = 0, is the fraction $ of the length over which the polarization 
r] takes a certain sign. Its PDF is p$((/>) = [ir^/ 0(1 — 0)]""^ [T7] . 

For b = {a{x)) 7^ 0, if 6 has opposite sign to the ellipticity angle 

T]o = sin-i [tanh (ln(|c/+|/|rf_|))] /2 

of the injected pulse, the soliton polarization will first become linear after a distance X^^, 
move near the favored circular polarization over a subsequent distance X^ra., and never leave 
this ultimate circular polarization after a subsequent distance Xq^c- Thus, X^^p = + 
Xtra + Xfluc. The first two distances are distributed as |T7] 

(8) pxix) = (^\b\i/aV2'Kx^^ exp [-b'^{£ - xf /2a^x\ , 

and (X) = £, (t| = a^f/fe^, with 

l = \n{\d+\/\d.\)/2\rb\ 



^(4n + 2)e" 



(4n+2)2L/2a; 



16n2L/2x 



n=0 



n=l 



for X = Xlin and 



£ = tanh"^(cos2K)/2|r6| 
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for X = Xtra- The probability density of Xfjuc is 



(9) Px^Jx) = {\b\/aV2^) exp{-b^x/2a^). 



When b and r/o have the same sign, results are similar. The above two distributions show 
excellent agreement with the numerical simulations in the right panel of Fig. |3l 

Numerical simulations of system ([1]) are carried out by aligning a grid in x and t with 
spacing Ax and At to the characteristics of Eq. f|Ta|) (i.e. Aa; = At) and implementing a 
scheme based on the implicit mid-point method. The initial conditions in n± are obtained 
using Eq. (4) in which a{x) is a Gaussian random variable at each grid point with mean b 
and standard deviation a/y/Ax. The soliton enters through the boundary condition on the 
left hand side, and no boundary condition is needed on the right hand side. 

In conclusion, we have shown excellent comparison between the analytical results and 
direct numerical simulations describing the statistics of the orientation, ip, and eccentricity, 
7], of the light polarization ellipse. We also found favorable comparison for the distance 
over which the light polarization evolves from linear to nearly-circular, Xtra, and back again 
to linear, Xjnt, and, for unequal distributions within the lower populations, the distance 
between the point at which light first comes close to its favored circular polarization and the 
point at which it finally forever remains near it, Xfluc- 
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